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1 Introduction 

Over the past several decades there have been remarkable breakthroughs 
in experimental techniques for probing non-classical behaviors of small-scale 
quantum objects coupled to quantum environments (see, e.g., [5]). Corre- 
spondingly, a more sophisticated theoretical understanding of the thermody- 
namic nature of such systems, especially in the low temperature regime where 
quantum effects are dominant, has been substantially demanded. In contrast 
to common quantum statistical mechanics, which is intrinsically based on a 
vanishingly small coupling between system and environment, the finite cou- 
pling strength between them in the quantum regime causes some subtleties 
that must be recognized. 

At the heart of the aforesaid "quantum thermodynamics" [5j i the second 
law of thermodynamics, assumed to be inviolable by the scientific community 
for over a century jl], has been confronted by challenges with considerable 
interest, and its absolute status has even come under increased scrutiny 0I3J 
[5]. In fact, this fundamental law of nature has extensively been theoretically 
studied particularly in the scheme of a quantum harmonic oscillator linearly 
coupled to an independent-oscillator model of a heat bath (quantum Brow- 
nian oscillator) in equilibrium at a (low) temperature T, mainly due to its 
mathematical manageability. 

A short overview of recent results either pro or contra the validity of 
the quantum second law has been given in [T|; and the final result therein 
was that a properly defined Clausius inequality 6Q* S < T* s cISn represent- 
ing the second law is valid without any previously argued violation in form 
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of SQ ^ T dS at T — > 0. In the above inequality, represents the von- 
Neumann entropy of the coupled oscillator, which is, interestingly enough, 
identical to the thermal entropy of an uncoupled effective oscillator in equi- 
librium. And SQ* S is a heat exchanged between (weakly coupled) effective 
oscillator and bath, and T* ff a well-defined effective equilibrium temperature. 
This effective temperature differs from the temperature T of the total system 
(oscillator plus bath) especially in the strong-coupling limit, where the total- 
system temperature is in fact not well-defined as an equilibrium temperature 
of the system oscillator since the reduced equilibrium density operator of 
the oscillator H s is not any longer in form of the canonical thermal state 
P/3 oc e~P Hs with (3 = l/(ksT) [T|. This discrepancy between these two 
equilibrium temperatures is, of course, caused by the extra work (or energy) 
needed to couple an (uncoupled) oscillator to a bath [B]. As a result, it may 
be legitimate to say that the original form of the Clausius inequality for the 
coupled oscillator in terms of the temperature T is not well-defined rather 
than being violated. 

However, the entire discussion of the quantum second law has so far been 
restricted to that for thermal equilibrium states, accordingly, in form of the 
Clausius equality. On the other hand, there has recently been an interest- 
ing result for a generalized Clausius inequality for non-equilibrium quantum 
processes, but restricted to isolated quantum systems (initially prepared at a 
thermal equilibrium state) [9] . In this paper we extend the above discussions 
into far-from-equilibrium processes in the scheme of quantum Brownian os- 
cillator as a prototype of open quantum systems, but in the weak-coupling 
limit only since the exact treatment of the non-equilibrium processes in the 
strong-coupling limit is pretty much hopeless to leading to useful quantities 



in closed form to be needed for our discussion [cf. Eqs. ( 28 )-( 36 1] . To do 
so, we consider the quantum oscillator with time-dependent Hamiltonian pa- 
rameters, which finally enables us to derive a generalized Clausius inequality 
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without any violation in variation of the Hamiltonian parameters as our cen- 
tral result. In fact, the time-dependent quantum oscillator has been studied 
by many authors for last about 60 years, initiated by Husimi 1 10] . In this 
paper we appeal to the method of quantum Liouville equation in order to 
directly obtain the time-dependent density operator of the system in consid- 
eration. This approach differs from that developed by Husimi, which is based 
on the fact that the Schrodinger equation for an isolated linear oscillator, 
rather than coupled to a bath, can reduce to a system of classical equations 
for the oscillator. 

The general layout of this paper is the following. In Sect. [2] we review 
the general results of quantum Brownian oscillator needed for our later dis- 
cussions. In Sect. [3] we consider non-equilibrium processes and then derive 
a closed expression for the time-dependent reduced density operator of the 
oscillator weakly coupled to a bath at an arbitrary time along the processes. 
In Sect. [4] the same discussion will take place for the corresponding quasi- 
static processes. Next the second law of thermodynamics for non-equilibrium 
processes will be systematically discussed in Sect.[5j Finally we give the con- 
cluding remarks of this paper in Sect. [6j 



2 Basics of quantum Brownian oscillator 



The quantum Brownian oscillator under investigation is described by the 
model Hamiltonian (Caldeira-Leggett model) [TT|[T2] 



Hq — H s + Hb + H s f, , 



(1) 
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where a system linear oscillator, a bath, and a system-bath interaction are 



i=i 



iV W 2 



respectively. Here the coupling strengths Cj, and the spring constants k Q 



Mloq and fc, = rrijUn. The total system _ffo is assumed to be within the canon- 



ical thermal equilibrium state pp — e~^ H ° /Zp, in form of a non-separable 
state (9^ p s £g) /)(,) due to the interaction H s b. where the partition function 
Zp = Tre^^ H °. The second term of the interaction H s b, proportional to 
q 2 , was introduced in order to protect the pre-determined frequency (wo) of 
the system oscillator H s from its modification induced by the system-bath 
coupling (the first term linear in Xj) Here the system and the bath effec- 
tively share the energy in the coupling term, especially in the strong-coupling 
limit (a,- ^> 1), and so it is in fact not completely clear whether this energy 
should be interpreted as belonging to the system or to the bath [T3]. There- 
fore, without the above second term, the internal energy U s = Tr (H s pp) 
of the coupled oscillator alone, with its unique frequency w , as well as its 



reduced density operator would not be well-defined [cf. Eqs. (59 1 and ( 65 1 ] . 
In fact, from the Heisenberg equations of motion for q and p we can derive 
the quantum Langevin equation without the frequency shift as |11[I12| 

Mq\t) + M f dT-y(t -r)4(r) + Mcu 2 q(t) = £(t) , (4) 
Jo 

where the damping kernel and the noise operator are, respectively, given by 
1 N c 2 

3=1 3 * 

i(t) = -Mf(t)q(0) + J2 C J co8(w,-t) + sin^t) j . (6) 
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Here as required, (£(t)} Pbl = Tr py} = 0, in which the shifted bath state 
P&/ = (Hb+H sb ) J ) ^^jj t ne corresponding partition function \ and 
the noise correlation [H] 

, N 2 

- &o + 1(0 !(*)) = ~ E — C0S K(* ~ 0} coth 

(7) 




Now we introduce a response function |12| 

XA B {t) = \ ([A(t),B(0)}) p O(t), (8) 

where 0(t) represents a step function. Then it can easily be shown that 
X Pq (t) = -Xqp(t) = Mxgq(t) and Xpp(t) = ~M 2 Xqq(t)- For a later purpose 
it is also necessary to discuss the time-reversal dynamics of q(t) in terms of 
f(t) := q(—t) and its momentum s(t) :— —p(—t). We can then derive the 
corresponding quantum Langevin equation |14j 

Mf(t) + M [ dTj(t-T)f(r) + Mul?{i) = (9) 
Jo 

While this is the same in form as Eq. Q. the two equations differ in the 
noise term in such a way that is identical to £(t), however, with replace- 
ment of pj (0) — > —pj (0) in ([6| . From Eq. Q and the stationarity relation 
(A(t)B(0)) p = (A(0)B(-t))p [12_ it appears as well that Xrr(t) = ~Xqq{t) 
and Xrs(t) = -Xrp(t) = ~Xgp(*) and Xss(t) = ~Xpp(*)- Applying the Laplace 
transform technique to Eqs. Q and ([9]), respectively, we can finally obtain 
the exact expressions [TS] 

«K*) = -X9p(*)e + Xgg(*)p-E^«w(*)^ (10a) 

i 

^(*) = -Xrp{t) Q + Xrr(t) P - *^2{Xr P] (t) X-j ~ Xr Xj {t) Pj} , ( l ^b) 

3 
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where the operators g,p, £j, and pj represent the initial values q(0),p(0), %(0), 
and pj (0) , respectively. Here we have 

X gXj (t) = ^ J djjx qxj {w)e-^ (11) 

with Xqxj (w) = Cj/{rrij (w| — w2 )} ■ X<?g( w )i where the susceptibility x gg (w) 
is the the Fourier-Laplace transform of Also, Xrxj(t) — —Xqx j (t)i 

Xrpj (t) = X 9P3 (t), and x gPj (*) = -m 3 - Xqx } (*)• 



It will also be useful later to introduce the well-known expressions for the 
equilibrium fluctuations in terms of the susceptibility x qq (uj) such as |16| 

= ^J duJ Coth (2) Im fe^ + l0+ )> ( 12a ) 

<p 2 )^ ~/" Wcoth (^r) Mx TC (w+<o+)}, (m) 

which can be derived from the fluctuation-dissipation theorem [17] . For the 
Drude model (with a cut-off frequency ujd and a damping parameter j a ), 
which is a prototype for physically realistic damping, the equilibrium fluctu- 
ations are explicitly given by |18| 

^iX>2'{4 + M^)} (13) 

<f)? = -Mp?^{^ + ^(^)}, (14, 

respectively, where the digamma function ip(y) = d In r{y)/dy [19|, and wi = 
J?, 0^2 = ^1,1^3 = 22, and the coefficients 

,(i) + z 2 .(2) __ O + z 2 (3) i7 + zi 



(12 - zi)(z 2 - n) 1 d ( Zl - Q)(z 2 - Zl ) ' d (z 2 - fl)(zi - 

(15) 
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Here we have adopted, in place of (w , 7o)j the parameters (w ,i2, 7 ) 
through the relations [2U] 

"° := wo 2 ^ ; c d := /7 + 7 ; 7o := 7 ^ j9_ , (16) 



and then zi = 7/2 + iwi and Z2 = 7/2 — iwi with wi = y (wo) 2 — ( 7 /2) 2 . 



From Eq. ( 16 ) it also follows that 

^ 3 -w rf r2 2 + (cj 2 + 7o w d )/2-w d u; 2 = 0, (17) 

which will be used later. And we can then obtain for this damping model the 
response function expressed as |15| 

XwU " M («-«i)(»i-«a)(«a-fl) ' 1 1 

which is real- valued and vanishes at t = and 00. It is also interesting to 
note that this response function is temperature-independent indeed, which 
was originally defined in Q as a function of temperature. In fact, the sus- 
ceptibility, defined as the Fourier-Laplace transform of the response func- 
tion Xqq(t)> is explicitly given by the temperature-independent expression 
X qq {oj) = — l/{M(cj 2 + iu; 7 (aj)— w 2 )} in terms of j(u>) defined as the Fourier- 
Laplace transform of the damping kernel 7 (i), which can easily be obtained 
by applying the Laplace transform to Eq. Q |12| . 



3 Non-equilibrium process and its reduced density operator of the 
coupled oscillator 

Now we disturb the system of interest by varying its Hamiltonian parameter 
with time, namely, either the spring constant k(t) of the coupled oscillator or 
its mass M(t). Therefore, we should deal with a time-dependent total system 
H(t) — H s (t) + Hi, + H s b from now on, where the time-dependent coupled 
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oscillator H s (t) is explicitly given by either 



~ + ^?=H s + h 1 (t)?=: 1 H s (t) (19) 



l '~ k °q 2 = H s + h 2 (t)p 2 =: 2 H s (t). (20) 



2M(t) 2 

Here the initial values k(0) — ko and M(Q) = M, and h\{t) — {k(t) — 
k }/2 and h 2 (t) = {M - M(t)}/{2M ■ M(t)}. To derive the (time-dependent) 
reduced density operator of the oscillator % s (t) , we first consider the equation 
of motion for the density operator of the total system, which explicitly reads 

as nuna 

p(t) = e- itLo p(0) - % [ dre-^-^^L^pir). (21) 
Jo 

For a variation of the spring constant, the total Hamiltonian is Hi(t) = 
Hq + hi(t)q 2 , and the corresponding Liouville operator = Lq + 
satisfies 

Uh(t) = \[H^Px{t)] ; L { 1\t) Pi {t) = l -[q 2 1 pi{T)]h 1 {T). (22) 



Here the Liouvillian L^\t) surely corresponds to L\(t) in (211. Likewise, 



for a variation of the mass the total Hamiltonian is H 2 (t) = Ho + h 2 (t)p 2 , 
and accordingly — L + with <— L\ and 

Loh® = \[H ,f>2(t)] ; L {2 \ T )p 2 (r) ^ l -[f,p 2 (T)]h 2 (r). (23) 



Now we attempt to obtain the density operator pi (t) in its explicit form. 
To this end, we mimic the technique applied for the study of field-induced 
dynamics in the quantum Brownian oscillator, discussed in |15| : we first sub- 



stitute (22 1 into (21), with pi(0) = pp, and then make iterations for pi(r) in 
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the integral. Then we can arrive at the expression 



+ {-{f Jldrh^T) dr' hx{T')e-^-^ \q 2 , e -^- T ')"o [q 2 ,pp] e ^ T '^« 
+ (-if Jo dT Mr) Jo M Mr') f T ' dr" h x {r") e~W-)^ 



,i(t-r)H 



e -%(T-r')H 



, e -i(r'-r")ffo e i(r'-r")ffo 



,i(r-T')H 



e t(t-r)H + ... .(24) 



With the aid of [/5,g,7?o] = 0, this equation easily reduces to the expression 
in terms of f{t) = e'Ti tH °qeTi tH ° as 

Pi(t) = h + mfo dTh i(r) [r 2 (t-r),pp] + 

(m) 3 Jo & Mr) Jo dr' h^r') dr" h x {r") [r 2 (t r), [f 2 (t r'), [r 2 (t r"), ft,]]] 
where [f 2 {t),pp] = f(t) [r(t),pp] + [r{t),pp] f(t) with 

[r(t),Pp] = -Xrp(t) [q,Pp] + Xrr{t)[f>,Pp]~^X rP] {t)[x j Jp]+^Xrx ] {t)[Pj,h] , 

3 3 

(26) 



(25) 



obtained directly from Eq. (10b I. Likewise, we plug (23 I into (21 1 and then 



apply the same technique as that used for (24 1, finally leading to the density 



operator p2[t) in its explicit form, identical to Eq. ( |25[ ) but with replacement 
of all hi(t) -> h 2 (t) and all f 2 {t) ->• s 2 (t) where s(t) = -e~% tA °peTi tAo . 



Also, from (26 1 and s(t) = Mr(t) we have 



K*)>P/3] = -Xs P {t) [q, pp] + Xsr{t)\p,pp] -^Xm W)[Xj, Ppl + ^Xsxjit) [p. 

j 3 

(27) 



It is also instructive to rewrite Eq. (251 as its compact form 



Pl (t) = Tcxpj^^ dThi(T)f 2 (t-T)\ pp exp j-^^ dTh!(T)r 2 (t-T)^ , 

(28) 
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which is equivalent to 



{-I 

\ih Jo 



Texp<j — / dr/n(r)P(t-r)| Up I Texpi-~J dr hi(r) r* (t - r 

(29) 

The ordinary time ordering operator T and its reverse time ordering T were 
introduced; in fact, f 2 (ti) and ? 2 {t2) do not commute at different times t\ ^ 



ti [cf. Eqs. (31al-(34cl]. Here we used the well-known operator identity 



XB Ae -XB 



= A + A 



B,A 



21 



B, 



B, A 



3! 



B. 



B. 



B, A 



(30) 



Let f 2 (t) = Si(t) + Bi(t) + Ci(t), where explicitly 



Si(t) ■= {Xr P (t) q - Xrr(t)p} 2 (31a) 

&l(t) ■^^{Xrp ] {t)x j -Xrx 3 {t)Pj}{Xrp k {t)x k -Xrx k it)Pk} (31b) 



Ci(i) :=2{x n >(t)g-Xrr(*)^}5^{Xrp^(*)*j-Xr!B i (*)Pj} , (31c) 



and the three operators surely commute pairwise at any time t. On the other 
hand, we have at ti =/= t 2 



5i(*i),Bi(t 2 ) = 



(32) 



but all other commutators of the operators do not vanish indeed, which are, 
respectively, in form of 

SiCii),^)] = J2 {fj^itlM qx 3 + ff\h,t 2 ) QPj + f?\t 1 M)px ] + f^(h,t 2 )ppj} 

3 

(33a) 

Bi (ti), di (t 2 )] - ^ {. 9j (1) (ti, t 2 ) qx s + gf (h , t 2 ) qpj + gf (ti , ta)^ + ffj 4) *a) #j } 
i 

(33b) 



12 



as well as 



<Si(*i),Si(i 2 ) =ai(ti,t 2 )q + a 2 (h,t 2 )p + a 3 (*i, i 2 ) {qp + pq} 

(34a) 

Bi(ti),Bi(t 2 )] = J2 {^fkih^XjXk + X^ih^XjPk + \${t x ,t 2 )PjXk + \?^(ti,t2)fijPk} 

(34b) 

£l{ti)A(t 2 j\ = *2) <f + Utl,t 2 )f + &(ti, t 2 ) qp + &(«!, t 2 ) pq + 

(34c) 



,,.fc 



Here /™ fa , * 2 ) , 5 n) (t x , t 2 ) , a„ fa , t 2 ) , A % ] fa , t 2 ) , £ n (t i , h) , £ S ; (ii , h) 



are some 



scalar functions. 



These non-commuting properties make it highly complicated to explic- 
itly carry out the transformation of the time-ordered exponential operator 



Te i/-r (t-r) in Eq ^ with the aid of Eqs . (32|-(34c) and (46l-(48| as 



well as the Zassenhaus formula and its dual, the Baker-Campbell-Hausdorff 
formula [22 , 23 , into its factorized form of 



(35) 



where (• • • ) = a,j(t)qXj + bj(t)qpj + Cj(t)pXj + dj(t)ppj in terms of the 
system-bath couplings only. In fact, this transformation process is a critical 
step for obtaining the reduced density operator lZi(t) := Tri, pi(t) of the 
coupled oscillator {H s (t) in its closed form from the total density operator 
Pi{t) in such a way that, by the cyclic invariance of the trace, 



Hi (t) = Tr b | e •^i=»"»-f-'e ■'-^ >e " J ^"->"' pp e 
o-i'EA- ) e ~ifm,t) * Jf(q,P,t) JT,A- ) 



(36) 



Trb {, <>J ) e -*/(M,t) ^ e i/(M,t) e »E,(- ) } 
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Here, Trj, denotes the partial trace for the bath alone. And the initial equilib- 
rium state 71(0) := Trt, pp is defined as the reduced operator of the canonical 
state pp and explicitly given by 



( q \n(o)W) 



(q + q 1 ? _ (P 2 )p (q - q'f 
%{q 2 )p 



exp 



2h? 



(37) 



which holds true regardless of the system-bath coupling strengths. Likewise 
£a(t) :=Tr fc p 2 (t). 

Consequently we now restrict our discussion for a closed form of the re- 
duced density operator !Z\(t) to the weak-coupling limit, where Xrxj (t), Xrpj (t) 



and so especially d{t) 0. From Eqs. (|28J, (|31c| and (|32|-(|36), it then 
follows that 



li[ w \t) wTexpj^ J dr/ii(r)<Si(f-r)| ^(0) exp|-i^ dr hi(r) &(* - r) 

(38) 

We stress here that this weak-coupling limit obviously differs from an isolated 
system with identically vanishing coupling strengths (cj = 0); in fact, the re- 
sponse functions Xrp(t) and Xrr(t) of <Si(t) depend on the coupling strengths 
already, as was discussed in Sect. [2] Also, it is worthwhile to point out that 



Eq. (38) can be regarded, by construction, as a good short-time approxi- 



mation to an exact expression of the reduced density operator lZi(t). Fur- 
ther, as the response function Xqq(t) m (18) and so the resulting quantities 



{Xr% (t), Xrp 3 ■(£)} [cf. exponentially decay with time, the contribution 

of C\ (t) to the density operator H,\ (t) may not be significantly non-negligible 
with time t large enough even in the strong-coupling limit, unless h\{t) ex- 
ponentially increases. 



Let us now simplify the formal expression of the reduced density operator 
(t) in ( 38 1 by considering, with the aid of ( 30 ) , its expanded form such 
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as (25 1; using Eq. (31a I we can first obtain 



(q\T n S!(t),Pf) W) = Ai(t)(q\K(0)\q'), 



(39) 



where Ai (t) := xL(t) {g 2 - (</) 2 }+2^ Xrr (<) Xr P (t) {qd q +d q . q')+{ih) 2 x 2 rr (*) (d 2 



Likewise, from (27) and s 2 (i) = S 2 {t) + B 2 (t) + C 2 (t) with S 2 (t) 
{Xsp(t) q — Xsr(t)p} 1 we can a l so have 



{q\Tx b S 2 (t),pp W) = A 2 (t) (q\TZ(Q)\q') , 



(40) 



in which A 2 {t) := X %{t) {q 2 -{q') 2 }+^h X sr(t) Xs P {t) {qd q +d q ,q , )+{iH) 2 X 2 sr (t) (d 2 
d 2 ,). Here we used Tr b ([xj,pp]) = Tr b Qpj,pp]) = 0. Similarly, Tr b (£j \p, pp]) = 
Tr 6 (pj [q, pp]) = Tr b (pj [p, p^]) = 0. And 

( 9 | Tr b {g, [g, \q') = {q 2 (q') 2 } (q\£(P)W) 

(q\ Tr b ({g, % pp]} + + {p, [q, pp]} + ) \q') = -2th (qd q + d q ,q') (q\K(0)\q') 

(q\ Tr b {p, \p,pp]} + \q>) = (ihf (d 2 - d 2 ,) <g|£(0)|g') , (41) 



where the anticommutator {A, B} + = AB + BA. 



With the aid of Eqs. ( 39 )-( 41 1, we can next obtain 



(g|Tr b [&(*), [^(rj.ppjj \q') = ii(i)ii(r) {q\K(0)\q') (42a) 
(g|Tr b [s 2 (t), [S 2 (r),p ]] W) = A 2 (t)A 2 (r) (q\n(0)\q>) . (42b) 

Along the same line, after making a lengthy calculation, we can also arrive 
at the expressions 

(q\Tr b [&(t), [&(r), [^(t'),^]]] \q') = A x {t) A x {t) A 1 {t') <g|£(0)|g') (43a) 
(g|Tr b \s 2 (t), \s 2 (r), \s 2 (t'), pp] ] ] \q') = A 2 (t) A 2 (r) A 2 (r') (q\K(0)W) -(43b) 
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Based on Eqs. (25 1, ( 39 1-( 40 1 , and (42a I- (43b I we can finally find the matrix 



elements of the reduced density operator of the coupled oscillator {H s (t) as 



(q\TZ^\t)W) = Te-i J ^) {q \TZ(Q)\ q ') 



(44) 



where J\{t) = J dr h\ (r) A.\ (t — r) represents the time-evolution action. 
Likewise, the density operator of the coupled oscillator 2 % s {t) can be obtained 



(q\H 2 w) (t)\q') = Te-^{q\1l(0)\q') 
where J 2 (t) = J Q dr h 2 (r) A 2 (t - r). 



(45) 



Now we see from Eq. (34a I that Ai(t\) and A\{t 2 ) at t\ ^ t 2 are not 



commuting and accordingly it is non-trivial to directly deal with the time- 



ordered exponential operator in (44). Therefore, we need to introduce the 



exponential operator identity derived in 



Texpj^ dr£ t (r)j = expj^# n (t)j , 



(46) 



where Bt(r) :— h\(r) A\{t — r), and the low-order terms are explicitly given 
by 



K 1 (t) = C 1 (t) ; K 2 (t) = ~C 2 (t) 
K 3 (t) = ~£s(t) + ^ [c 2 (t),C 1 (t)] ; K A (t) = \CS) + - [c^tlC x {t) 



(47) 



Here the commutators 



C n (t) 



drx dr 2 - ■ ■ dr n 
Jo Jo Jo 



B t (n), [B t (r 2 ), , [B t (r n _i),B t (r n ) 

(48) 

with C\(t) — J *drS t (r). In fact, the operators K n (t) for all n can be eval- 
uated exactly. 
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Let us simplify the commutators C n (t) to derive the closed expression of 
(q\TZ[ w) {t)\q') in ([ill. First let y:=q + q' and z := q - q', and so d q + d q > = 



2d y and d q — d q / = 2d z . Then it easily appears that Ai(i) = {Xrp(t)y + 
2ihxrr(t) d z }{xrp(t) z + 2ihxrr(t) d y } . This allows us to finally obtain 

Cn(t) = yz + b n (t) (yd y + d z z) + thc n (t) d y d z , (49) 

where a n (t) , b n (t) and c n (t) G K, and explicitly given by 

Q-2m-l(t) — 2 3 " 1 3 l2m-l,l(t,t) ; CL2m(t) = 2 3 ™ 1 -f2m,l(*)*) 
&2m-l(*) = 2 3 ™ 2 l2m-l.l{^^) > fem(*) = 2 3 " 1 1 ^2m,2(^) 

C2m-iW = 2 3m - 1 J 2m _ 1)3 (M) ; c 2m (t) = 2 3m+1 / 2m , 3 (t,t) (50) 

with m = 1, 2, 3, • • • . Here 

:= dr hi{r) xlpit - t) ; I 12 (t,u) := dr ht(r) Xrr(t - r) x rj) (t - r) 
J 13 (t,u) := r u dT^(r) X 2 r (i_ T ). (51) 

For n > 2, we have 



(52) 



and 



h m ,l(t,u) := J™dThi(T) {Xrr(t ~ t) Xrp(t ~ T") 7 2 m-l,l(*! T ) ~ Xrp(* - 1~) I 2m -l,2(t,T)} 
hm,a(ti u ) := fa drh l{ T ) {Xrr(t - T~) hrn-lAt^) ~ Xrr(t - t) Xrp(t - t) hm-l^' T )} 

(53) 
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as well as 

hm+l,l(t, u) 
^2m+l,3(*i u ) 



J Q U dr fti(r) {2 Xrr(t - r) Xr P (i - t) I 2m ,i{t, r) - Xr P (* ~ r ) hmA^ r )l 

Jo" G?T ^l(r) {Xrr(* ~ T ) h-mA*' T ) ~ 2 Xrr(t - t) Xrp(t - t) hm,z(t, t)} . 

(54) 



Here we employed the commutators in ( 58 1 with the replacement of M 



yd y + d z z and A — > yz and B — > d y d z , which immediately leads to a — > 
2, to — > — 1 and r — > 0. With the aid of Eqs. (46l-(49) we can then rewrite 



the time-evolution in (44 1 as the unitary operator 
T exp < — 



H ■*■«>} — {"f 



where the real-valued coefficients 



yz + b\ (i) (ydy + d z z) + ih~ c\ (t) d y d. 



(55) 



5i(t) := ai(t) + ^ + ^ + I {ai(t) b 2 [t) - h(t) a 2 (t)} + 
^l + l{a 1 (t)b 3 (t)-b 1 (t)a 3 (t)} + --- 

h(t) := h(t) + ^ + ^ + ~ {oxWoaW - Cl (t)a 2 (t)} + 
^ + ^ {ai(t) Cs(t) - ci(t) a 3 (t)} + • • • 

**(*) := ci(i) + ^ + ^ - I {d(t) 6a(t) - ^(t)} + 



C4(t) 

4 



{d(t) 6a(t) - 6 X (f)cs(t)}- 



(56a) 



(56b) 



(56c) 



(in fact, all higher-order terms can be determined exactly). Likewise, coef- 



ficients 02(i),&2(^) and c 2 (£) pertaining to Eq. (45) can also be introduced, 



which are identical to their counterparts in (56al-(56cl, respectively, however 



obtained from the replacement of hi — > /i 2 and Xrr ~^ Xsr & n d Xrp ^ Xsp m 
d5B-(|54h. 
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To further proceed with ( 55 1 , we apply another operator identity, derived 
in [23], given by 

exp | A (M + fiA + 2/S) j = cxp (kM) exp ^fe~ aK Aj exp (gBj exp (d) 

(57) 

where A, /i, v are arbitrary complex numbers, and the product \xv is assumed 
to be real-valued. Here the operators M, A, B satisfy the commutator rela- 
tions 



A,M 



aA ■ B,M = aB ; A, B = mM + r (58) 



where a, m, r £ R. And numerical functions / = \iX and g = vX where 
X = {ttm\D)/{D - (a/2) tan AL>} with D 2 = -{pvm + a/2)(a/2), and 
« = -(2/a)ln{cosAD - (a/2£>) sin AD} and d = (« - A)r/m. 



After making a lengthy calculation with the aid of Eq. ( 57 ) , every single 



step of which is provided in detail in Appendix, we can finally arrive at the 
closed expression 



1 gift) 



cxp -Bi{t) 



I 8<<f), 



) 2 , (tf)p{q-q'Y 



9 fl 

q z — q z 



(59) 



where the two dimensionless parameters 

D(t) 



B 1 (t) := 



D(t) cos{6i(£) D(t)} - sin{6 1 (t) D(t)} 



*i(*) :=-^fl(*)Bi(t) -%(?hf(t) e 



(60) 
(61) 



in terms of the coefficients a\(t),bi(t) and ci(i) in (56a)-(56c). Here, the 
parameter D(t) = ±{5i(t) ci(i)/6 2 (t) - 1} 1/2 el or iM as given in Appendix 
and so f(t) and g(i) s £R in ([82]) . As shown, the time-dependency of the 



reduced density operator in (59) consists entirely in B\{t) and ^i(^). Figs. 
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[l|2] demonstrate their behaviors versus time for, e.g., k(t) = k (2 — e _t ) 
within the Drude damping model. Applying exactly the same technique, we 
can also derive the reduced density operator (q\1Z^"\t)\q') in closed form, 
which is in fact identical to Eq. ( [59] ) but with replacement of B\ (t) , $1 (t) — > 
B 2 {t) , ( I ) 2(t) in terms of 02(^)1 and Ci{t) for f(t),g(i) and D{t) therein. 
The normalization Tr 7Z[ W ^ (t) — Tr7l%\t) = 1 can easily be verified with 
the aid of the identity [T5] 



&e-'" 2+2ta ) = J-eV. (62) 



a 

It then follows that (q) Pl ( t ) = (g) n M^ = and (p) Pl ( t ) = iP) n M(t) = °' 
and 

/ 9> 1 ;.9i /„9v / n , ~9i 4/l 2 ^?(i) 1 

(63) 

From this, the instantaneous uncertainty relation also follows as 

= (fh (f)p + ■ (64) 

Then the instantaneous internal energy of the coupled oscillator reads as 

i*4(t) = {<H s {t)) <){t) = <^g$a + *f <<f } Pl (<) . (65) 
Along the same line, the expectation values for the density operator TZ^ (t) 



easily appear, respectively, as the counterparts to those in Eqs. (63 1 and 



(64 1 in terms of B2(t) and @2{t), and so the instantaneous internal energy 
2lA s (t) — (2H s (f))„w, , will immediately follow as well. 

K 2 {t) 

Comments deserve here. The compact form of the density operator iZi"\t) 



in (59), valid for an arbitrary variation of the spring constant k(t), was 



clearly derived for the special initial condition 1Z(0) in (37 1, or equivalently, 



the canonical thermal equilibrium state pp of the total system Hq with 
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[pp,H a ] = 0. This then gave rise to the significant simplification in form 



in the step from (24 1 to (251, which subsequently led, with the useful re- 



lations in (41 1, to Eq. ( 44 1 and finally Eq. ( 59 1 . In the general case of the 
initial condition, on the other hand, it is mathematically not straightfor- 
ward to obtain an explicit form of the density operator TZ^"\t). Clearly, the 
time-dependent coefficients {a\(t), b±(t), ci(t)} in (56a)-(56c) are fundamen- 



tal ingredients to the time-evolution operator in (55) and so the reduced 
density operator W^[t). Then, as shown in (50l-(54|, the coefficients are 
expressed in terms of the parameter hi(t) representing the variation of the 
spring constant as well as the response functions Xrr(t) an d Xrp(t), defined as 
the average values with respect to the initial condition pp but temperature- 



independent indeed [cf. (18l] and reflecting the characteristics of the bath 
coupled to the oscillator in consideration. 

It may also be worthwhile to point out that Zerbe and Hanggi derived a 
master equation for the reduced density operator 7ti(t), however, restricted 
to i) the periodic potential, hi(t) q 2 —> (m/2) ■ e cos(i?i + ip) q 2 ; ii) the Ohimic 
damping; iii) the initial state of the total system given by an uncoupled one 
Po = Ps ® Pb EZ]) whereas this is obviously not the case in our study. Ac- 
cordingly, the initial state po cannot represent a thermal equilibrium of the 
coupled total system (oscillator plus bath), which is necessary for the dis- 
cussion of the Clausius inequality in Sect. [5] Further, in a damping model 
without cut-off frequency (such as the Ohmic), which is not physically real- 
istic, the validity of the second law in the quantum Brownian oscillator may 
not be guaranteed |18ll2"5] . 



4 Quasi-static process and its reduced density operator of the 
coupled oscillator 

For comparison with the above non-equilibrium processes, we discuss the 
corresponding quasi-static processes. Here the system of interest undergoes 
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change infinitely slowly and so remains in equilibrium exactly in form of Eq. 



(37 1 in every single step such that for any spring constant fc, 

(q + q') 2 (P 2 )0{k)-(q-q'f 



{q\R cq {k)\q') = 



1 



exp 



2h 2 



(66) 



(valid for an arbitrary system-bath coupling strength indeed) , where the ini- 
tial values (q 2 )g(ko) = (q 2 )g and (p 2 ) ,3 (fco) = (p 2 ) g. Apparently, this density 
matrix looks different from (q\R!\"\t)\q') in ( |59| ) , and in general not in form 
of a canonical thermal state oc e~P 1 ™ s ^ [T]. Eq. (59 1, however, reduces to 



its quasi-static counterpart in (66 1 indeed if k(t) — > at every single mo- 
ment: As demonstrated in Fig. [2] for the parameter <&i(t) of (59 1, where 
k(t) = fco (2 — e _i ) and so fc(oo) — > 0, we have <?i(oo) — > 0. From this and the 
initial value <Pi (0) = 0, it must follow that if k(t) remains infinitesimally small 
at every single moment, then &i(t) — » always. This immediately leads to 



(P 2 ) pi (t) B^t) (p 2 )g in (631, and then (<f) Pl(t) = (q^g/B^t) -> (q 2 )p(k) 



as well as {p 2 ) Pl (t) (p 2 )p{k) in (59). As a result, we can arrive at Eq. (66 1. 



Consequently, without any harm we can straightforwardly adopt here, 
with fco — > fc, all results for the initial equilibrium state i? cq (fco) = 7?.(0) 
obtained in [1.; we can introduce an uncoupled effective oscillator 



2M* ff (fe ) 



in the same state R eq (k ), with its internal energy U* s (k ) := (H* s (ko)) R eq (k ), 
being identical to the internal energy U s (ko) := (H s ) R eq (k ) of the coupled os- 
cillator H s , as well as its von-Neumann entropy S* s (ko) = —kg Tr{i? cq (fc ) lni? oq (fc )} 
k B {v(k Q ) + 1/2} ln{w(fc ) + 1/2} - k B {v{k ) - 1/2} ln{v(fc ) - 1/2} in terms 



of v(ko) = \/(q 2 ) p(ko) ■ (p 2 )p(ko)/h, being identical to that of the coupled 
oscillator. Here the mass of the effective oscillator is given by M* ff (fco) = 
(P 2 ) p ( k ) I U s ( fc ) , and the effective spring const ant k* e ( fc ) = ( fc + (p 2 ) p ( fco ) / { M ( q 2 ) g (fc ) } ) /2 . 
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Subsequently the effective frequency easily follows as 




-eVM = ^/S^4 + ^ (68) 



which also allows us to have 



f/* fl (fc ) = ^eff(fco) yj {q 2 )p{ka) ■ (p 2 )fj{ko) ■ (69) 

Therefore, for the single state R cq (ko) we now have two different pictures of 
the Hamiltonian in consideration, namely, the coupled oscillator H s (ko) and 
its uncoupled effective counterpart H* s (k ). 

Then it can be shown that the effective picture H* s (ko) is, remark- 
ably enough, exactly in the canonical thermal equilibrium state i? cq (fco) oc 
exp{-/3 c yfc ) • H* s {k )}, where /3* ff (fc ) = l/{k B T* s (k )} with the well- 
defined effective temperature r* ff (fco) = haj* s (ko)/(kB \n{l/^(ko)}). Here 
£p(k ) = {v(k ) - l/2}/{v(k ) + 1/2}. From this, it also follows that 

<? 2 >mm = 2 * 77 , coth( ^( fc o)-^: ff (fco) | 

^ = W ^ . ( 70b ) 



As a result, for the quasi-static process (66) we can take all expressions from 
Eq. ( |67| to ( |70b[ ) simply with replacement of k — » e.g., the internal en- 
ergy U* s {k(t)} = U s {k{t)} = (p 2 ) {k(t)}/2M + k(t) (q 2 ) p{k{t)}/2, which is 
surely different from its non-equilibrium counterpart \U S (t) in ( |65[ ) (note that 
the time-dependency of the quasi-static quantities comes entirely through 
the k- value specified by time i). Needless to say, in case that the coupling 
constants Cj — > 0, then H* s (k) — > H s {k) as well as T* s (k) — > T. Also, for 
the upcoming numerical analysis it is useful to point out that in the Drude 



damping model we substitute ojq — > u) (k) = k(t)/M into Eq. (17 1, which 



will give the expression of the parameter f2(k) in terms of {w(fc), 7 Q }, and 
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then those of j(k) = uj d - Q{k) and wg(Jfe) = {k(t)/M}{uj d / fl{k)}, respec- 
tively. And Zi(fc) = j(k)/2 + ixvi(k) and z%{k) — j(k)/2 — iwi(fc), where 
Wl (fe) = V{wo(fc)} 2 -{ 7 (fc)/2} 2 . 



It is also interesting to consider a temporal behavior of a distance between 
the non-equilibrium state TZ^"\t) and its quasi-static counterpart R eq {k(t)}. 
To do so, we adopt a well-defined measure Dl(t) := Tr (^j. (*) — R C q{k(t)}) 2 
[29j . which is, independent of the dimension of the Liouville space, between 
and 2. With the aid of 1191 



/- oo ^dg'exp{-a(g + g ') 2 -%-e / ) 2 + ^ 2 -9 /2 )} = (71) 



we can obtain 



Dljt) = !L I _ 1 +— 1 ) - (72) 

VW*) {q*)p{k(t)} + (<P)p) (P 2 )p + (P 2 ) {Ht)}) + 4fc 2 *?(t) <?%{*:(<)} ' 

In Fig. [3] this measure for fc(i) = fe (2 — e *) is demonstrated for different 
parameters. Similarly we can also have 

Mf4-\ -DCnM2 _. S 2Hy/B 1 (t) 



Di(t) :=Tr{Kr>(t)-K(0)} 



Dl(t) :=Tr{4 q {fc(<)}-i? cq (fco)} 2 = 



h I i i 

2 



(73) 



v / ((9 2 >/3 + (<fMM*)}) ((p 2 )^ + (p 2 MM*)» 

(cf. Fig. [4]) . Finally it should be stated that all results in Sect. [4] also hold for 
the density operator TZ^ (t) for the mass variation, simply by replacement 
of the subscripts 1 — > 2 and k(t) —> M(t) of all pertinent parameters. 
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5 The second law of thermodynamics 



Based on the results found in the previous sections, we will explicitly discuss 
the second law of thermodynamics in the quantum Brownian oscillator. To 
address this issue, we need first of all the first law of thermodynamics 

dU s = C®™ d P" + P» dE ^> ' ( 74 ) 

n 

where ^2p n dE n — r Tv(pd'H s ) — 5W S corresponds to an amount of work 
on the coupled oscillator, and ^2E n dp n = Tr("H s dp) = 5Q S an amount of 
heat added to the oscillator [3U] . Next we consider a specific non-equilibrium 
process (I), leading to a finite (and so experimentally measurable), rather 
than infinitesimal, change in those thermodynamic quantities, in which the 
system begins and ends in thermal equilibrium states but is driven away from 
thermal equilibrium at intermediate times. Then an amount of the work along 



the process starting with the initial state (37) is given, with no harm, by 



,1V.(/) = fdrk r-^ 1 ) = \ fdTk(e) pi{T) (75) 



cf. Eq. ( 63 1] . Note here that at the end point r = t, the system {H s (t) may 



not necessarily be in an equilibrium state but relax to the end equilibrium 



state R cq {k(t)} in (66 1. However, no work is performed during this final stage 
of thermal relaxation. Then the second law in its Kelvin-Planck form [4_ that 
this work cannot be less than its quasi-static counterpart is expressed as 

iW a (*) > W s {k(t)}, (76) 

where the work along the quasi-static process 



t 

2\ 



W s {k(t)} = f drk/ 9 ^!) = \ I ,/r/,{y-),.{/,(rl}. (77) 

J ° \ I R cq {k(r)} J ° 
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Fig. [5] demonstrates the validity of this inequality and so that of the sec- 
ond law. Notably, however, based on the fact that the equilibrium density 
operator R cq {k(r)} is in general not in form of a canonical thermal state 
for the coupled oscillator {H s (t) in consideration (rather than its uncoupled 
counterpart H*g{k(r)}), it can easily be shown that the quasi-static work 
W s {k(t)} cannot be interpreted as a well-defined free energy change of the 
coupled oscillator {H s (t) where < r < t. 

Here it is also worthwhile to shortly point out that there is an alternative 
formulation based on the partition function Zi{k(r)} — Tr e - ' 3 ^h t ) /Tr^ e~^ Hb 
where f3 = T) and the total Hamiltonian "Hi(r) = {H s (t) + + H s i, 

P1|T5CT^|^|52P5P3] . This immediately leads to the well-defined free 
energy ^{^(t)} = — lnZi{fc(r)}//3. As discussed in detail in [T] (the last 
paragraph of Sect. 3 thereof), however, the free energy J-"i{£;(t)}, containing 
by definition the coupling-induced (H s i,) contribution, is not valid for the 
coupled oscillator {H s (t) alone. 

Next we discuss the second law in terms of heat. To do so, we first take 
into account the internal energy U s {k(T)} = U* s {k(T)} °f the coupled os- 
cillator {H s (t) as well as its uncoupled counterpart H* s {k(r)}. The first 
law of thermodynamics then tells us that the internal energy change along 
the quasi-static process is [/^{^(t)}^ = Q s {k(t)} + W s {k(t)} = Ql s {k(t)} + 
W* s {k(t)} , which is tantamount to iQ c s(t) + {yV c g(t) along the corresponding 
non-equilibrium process (I) above. Here the non-equilibrium effective work 
iW ff(t) and its quasi-static counterpart W* s {k(t)} can be obtained directly 



from Eqs. (75) and (77 1, respectively, with replacement of the coupled oscil- 



lator {H s {t) by its counterpart H* s {k{T)} such that 

lVM*) = \f Q drk(r) (V) pi(T) 9 * M *J k ( T)} + ^W) d k K s {k(r)}^ (78) 



effl 



W: s {k(t)} = i / drk{r) ( {p 2 )p{k{r)} • 9 * M * {k ( T)} + <? 2 M%)} ■ d k K B {k(r)} 
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where 



(f)p(k) k(t) d k (f) p (k) (f)p(k)-d k {p 2 }p(k) 



d k k* a {k) 




{e)p{k) (q%(k) 



[note the discussion just before Eq. ( |68| with replacement of k — > k\. And the 
quasi-static effective heat Q* s {k(t)} can be expressed as J* dr k(r) T* s {k(r)}- 
9 k Sj\r{k(T)} in terms of the well-defined effective equilibrium temperature. 
Here the von-Neumann entropy S]\[{k) is identified with the thermal entropy 
of the effective oscillator as 

sw*>-^ + ^ !«■>'<*> 

(80) 



^{l^(*,}H 



[cf. Eqs. (70a I and ( 70b l] 



Now let iWs-eff(i) := iW s (i) — iW e ff(t), which can be interpreted as the 
work needed for "switch of picture" from the uncoupled effective oscillator to 
its coupled counterpart along the non-equilibrium process (I), and its quasi- 
static counterpart W s . e s{k(t)} :— W s {k(t)} — W* s {k(t)}. Substituting these 



two work functions into Inequality (76 1 and applying the above first law, we 



can immediately derive a generalized Clausius inequality 

lQeff(f) < QtsiHt)} + A-effW , (81) 



where i/WO) := lWs-eff(i) - W s . cS {k{t)} with iZA s _ ef f(0) = 0. Therefore, in 
the picture of effective oscillator we hold the standard form of the Clausius 
inequality in terms of the well-defined (effective) temperature, but with the 
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additional term iA a _ c g(t). This inequality can be considered as a consistent 
generalization of the Clausius equality <5iQ c ff = T* s dkSN valid for the quasi- 
static process, introduced in [T|. Obviously, the extra term iA s „ e ff(t) iden- 
tically vanishes in this case. And in the vanishing coupling limit (cj — > 0), 
where T* s (k) — > T as well as both !>V s _eff(<) — > and W s . c s{k(t)} — > leading 
to iZ\ s _ c ff (t) — s- 0, we can easily recover the ordinary form of the Clausius in- 
equality in terms of the equilibrium temperature of the total system. In fact, 
in the high-temperature limit, where the thermal fluctuation in the coupled 
oscillator is predominant to the system-bath coupling H s i>, the ordinary Clau- 
sius inequality follows as expected. In the low-temperature limit, on the other 
hand, the additional term iZ\ s _ e ff(i) may not be neglected |35) ; in Fig. [6] we 
compare the incomplete Clausius inequality iQ c s(t) < Q* s {k(t)} (with no 



violation) with its complete counterpart in (81 ). As a result, we see that In- 



equality (81 1 is a generalized Clausius inequality representing the second law 
in the quantum Brownian oscillator, without any violation. Finally it should 
again be stated that all results in Sect. [5] also hold for the density operator 
72-2 (t)i simply by replacement of the subscripts 1 — > 2 and k(t) — > M(t) of 
all pertinent parameters. 



6 Conclusion 

In summary, we have analytically studied non-equilibrium dynamics in the 
quantum Brownian oscillator and then systematically discussed the second 
law of thermodynamics. We have first derived a closed expression for the 
time-dependent reduced density operator of the coupled oscillator in the 
weak-coupling limit along the non-equilibrium process. Based on this den- 
sity operator, we have found a generalized Clausius inequality in terms of 
the "effective" parameters, which indisputably reveals the robustness of the 
second law in the quantum regime. In introducing the effective picture, we 
reasonably required all thermodynamic variables to exist and to obey the 
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basic relationships, especially the first and the second laws. This method, 
as given, works for the harmonic oscillator but cannot easily generalized to 
apply to a broader class of quantum systems. Therefore the question about 
the (rigorous) validity of the second law for such systems remains open. 

Our finding can be considered as a consistent generalization of the Clau- 
sius equality valid for the quasi-static process, introduced in [T]. We believe 
that this inequality will provide a useful starting point for later useful discus- 
sions of quantum thermodynamics and quantum information theory within 
the quantum Browian oscillator as a prototype of quantum dissipative sys- 
tems; as an example, a consistent quantum generalization of the Landauer 
principle representing the computational irreversibility may be in immediate 
consideration, which can be understood as a simple logical consequence of the 
Clausius inequality [36 , 37 . Lastly, it is also desirable to numerically study 
the non-equilibrium dynamics in this system in the genuine strong-coupling 
limit as the next task. 
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A : Derivation of the density matrix in Eq. (59) 



From comparison of Eqs. (551 and (571, it easily follows that M — > yd y + d z z 
and A — > yz and B — > d y d z , and so a — > 2 and m — > — 1 and r — > 0. Let 
A — > b\(t) and /i — > 5i (t)/ih b\ (t) and v — > ihci(t)/bi(t). Then we have d — > and 
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D 2 -+ oi(t) ci{t)/bl{t) - 1 and 

1 5i(t) tan{foi(t)L»} 



/ 



»fi &i(t) D-tan{6i(i)D} 



g -H»f-# J^ffl^L (82) 



6i(t) D-tan{&i(t)D} 



k — > In 



D cos{6i (t) £>} - sin{fei (t) D} 



(83) 



With the aid of (57 1 and (581, Eq. (55 1 can then be rewritten as 



T exp \ — ^ J\(t) \ = exp(Kydy) exp(K,d z z) exp (/e 2K yz) exp(gd y d z ) . (84) 



Substituting (841 into (44 1 and then using the identity e cx " a: F(x) = F(e c x) for a 



function F |23| and d z z — 1 + zd z , we can easily obtain 



( q \7Z[ w, (t)\ q '} = (f e"») ■ {e a9 * d * (q\n(0)\ q ')Y^2 



(85) 



where the initial equilibrium density operator in (371 is then expressed in terms of 
y and z as 



(q\H(0)\q') 



exp 



1 2 {P 2 )l3 2 



2ft 2 



(86) 



Let e 9d y a * (q\n(0)\q') = e u9 i a * exp (-y 2 - ~z 2 ) =: G(y, 5) in Eq. g with the 
aid of g, where j/ := y/ ^8 (g 2 )^ and 5 := z^/{p 2 ) p /2h 2 , and u := ( 9 /4fi) ^ (p 2 ) p / (q 2 ) p ■ 
We subsequently consider the expansion 

5) = ^ — u n (dy d- z ) n exp (-y 2 - 2 

71 = 
OO ^ 

= ^ —j u" Hniy) H n (z) exp 



-2 -2 

y - z 



(87) 
(88) 



Here we used an identity of the Hermite polynomial, H n (x) = (—1)™ e x (d/dx) n e 
(I5j . The Mehler formula [381139] 



E— r tt n -ff n (a;i) H n (x2) 



s/l - 4u 2 



exp 



Aux\X2 — 4 (ai? + a: 2 .) w 2 
1 -4u 2 



(89) 
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then allows us to have 



From Eqs. (831, (851 and (90 1 we can finally obtain Eq. (591 
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Fig. 1: (Color online) -Bi(t) versus time t for spring constant k(t) = fco (2 — e ') 



and dimensionless temperature ksT/hwo [cf. Eq. ( 60 )] ; dashdot: T — 10 (high 
temperature), and solid: T — 0.01 (low temperature). Here h = kB = wo = M = 
ko = 70 = 1, and uid — 5. Since this spring constant exponentially decays with 
time, the coefficients 5i(t),&i(t) and ci(t) in |56al-(56cl should fast converge, and 



so the infinite sum over n in ( 46 1 can be replaced by a finite sum with the upper 
bound TV not necessarily large enough. This can be verified here by considering 
TV = 1,2,3 and 4. From top to bottom (at t = 10), (black dashdot: TV = 1), (blue 
dashdot: TV = 3), (red dashdot: N = 4), and (green dashdot: AT = 2); (green solid: 
N = 2), (red solid: N = 4), (blue solid: JV = 3), and (black solid: TV = 1). From 
this numerical result with an oscillating and fast converging behavior of Bi (t) with 
increasing iV, we adopt {-Bi(i)|jv=3 + Si(t)|jv=4}/2 as a numerical fitting of Bi(t), 
which will be used for later numerical studies. 



Fig. 2: (Color online) <I>i(t) versus time t for spring constant k(t) = fco (2- 



and dimensionless temperature ksT/hwo [cf. Eq. (61 1]; dashdot: T = 10 (high 
temperature), and solid: T = 0.01 (low temperature). Here % = kB = wo = M = 
fco = 7o = 1, and cj,j = 5. As for Fig.[l] we consider AT = 1, 2, 3 and 4. From top to 
bottom (at t = 10), (black solid: N = 1), (blue solid: AT = 3), (red solid: N = 4), 
and (green solid: A r = 2); (black dashdot: TV = 1), (blue dashdot: TV = 3), (red 
dashdot: TV = 4), and (green dashdot: TV = 2). From this numerical result with 
an oscillating and fast converging behavior of <I>i(t) with increasing TV, we adopt 
{^i(t)|iv=3 + $i(i)|iv=4}/2 as a numerical fitting of $i(t), which will be used for 
later numerical studies. Especially for T — 0.01, $i(t)|jv=3 ~ ^i(t)|jv=4 already. 



Fig. 3: (Color online) D%(t) versus time t for spring constant k(t) = fco (2 — e *) 



and dimensionless temperature fcsT/fiu;o [cf. Eq. (72|[; dashdot: T = 10 (high 
temperature), and solid: T — 0.01 (low temperature). Here h = kB = wo = M = 
fco =7o = 1- From top to bottom (at t — 10), (blue solid: cut-off frequency 
= 15), (red solid: cud = 5), and (black solid: Ud = 1); (black dashdot: cjd = 1), 
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(red dashdot: u>d = 5), and (blue dashdot: u d = 15). The three curves at T — 10 
are almost identical. 

Fig. 4: (Color online) {Dl(t), solid} and {D^.(t), dash} versus time t for spring 



constant k(t) = ko (2 — e ) and dimensionless temperature fcsT/fttJo [cf. Eq. ( 73 1] . 
Here h = kB = wo = M = ko = 7o = 1. Solid: from top to bottom (at t = 10), 
(green: T = 0.01 and oj d = 15), (black: T = 0.01 and u d = 5), (red: T = 10 and 
w,i = 5), and (blue: T = 10 and aid = 15). Dash: from top to bottom (at t = 10), 
(black: T = 0.01 and oj d = 5), (green: T = 0.01 and uj d = 15), and (red: T = 10 
and w<j = 5) « (blue: T = 10 and ojd = 15). 

Fig. 5: (Color online) iW s (t) — W^ s {fe(f)} =: j/5 > versus time t for spring 



constant k(t) = fco (2 — e ) and dimensionless temperature fesT/fiuo [cf. Eq. |76l]; 
dash: T = 10 (high temperature), and solid: T = 0.01 (low temperature). Here 
h = = (jJo = M = ko = 70 = 1- From top to bottom (at t = 10), (black dash: 
uid = 15), (red dash: w<j = 5), and (blue dash: ojd = 1) for 2/5/20, and then (black 
solid: uj d = 15), (red solid: aid — 5), and (blue solid: w<j = 1) for j/5. 

Fig. 6: (Color online) {I c ff (t) — Jes(t)}/ Jeff (*) — : J/6 versus time i for spring con- 
stant fc(t) = fco (2 — - e - ') and dimensionless temperature ksT/hoJo- Here T e ff(t) := 



Q* eS {k(t)} - iQeff(t), and J cff (t) := Q* ff {fc(t)} - ^(t) + iZi s _ eff (t) [cf. Eq. ([ST])], 
which is identical to iW s (t) — W B {k(t)} in Fig. [5} dash: T = 10 (high temperature), 
and solid: T = 0.01 (low temperature). And fi = ks = wo = M — ko = 70 = 1. 
Solid: from top to bottom, (blue: uj d = 15), (red: uj d = 5), and (black: ujd = 1)- 
Dash: in the same order. The three curves at T = 10 almost vanish, as expected. At 
T — 0.01, we see that i/A s _ e ff(i) < for uj d = 5 and 15, immediately leading to no 
violation of the incomplete inequality I e g(t) > 0; on the other hand, i/A s - e ff(i) > 
for Wd = 1. However, from the fact that 2/6 > —1 in this case, there is still no 
violation of I c ff(t) > 0. 
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